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Abstract: Peristaltic transport of a Bingham 
fluid in a channel with permeable walls is 
studied under long wavelength and low 
Reynolds number assumptions. This model 
can be applied to the blood flow in the sense 
that erythrocytes region and the plasma 
regions may be described as plug flow and 
non-plug flow regions. The effect of yields 
stress, Darcy number and slip parameter on the 
pumping characteristics are discussed through 
graphs. 
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I. INTRODUCTION 

Peristaltic transport is a form of fluid 
transport that occurs when a progressive wave of 
area contraction or expansion propagates along 
the length of a distensible channel/tube containing 
the fluid. Peristalsis is used by a living body to 
propel or mix the contents of the tube such as, 
transport of urine from the kidney through the 
ureter to the bladder, food through the digestive 
tract, bile from the gall bladder into the duodenum, 
movement of ovum in the fallopian tube etc. It is 
accepted that most of the physiological fluids 
behave like non-Newtonian fluids. This approach 
provides a satisfactory understanding of the 
peristaltic mechanism involved in small blood 
vessels, lymphatic vessels, intestine, ductus 
efferentus of the male reproductive tract and in 
transport of spermatozoa in the cervical canal. 

Nicoll and Webb (1946) and Nicoll reported 
that peristalsis plays an important role in blood 
circulation. Some theoretical and experimental 
investigations have been made on the peristaltic 
motion of blood considering blood as a non- 
Newtonian fluid. Latham (1966) first made an 
experiment to study the fluid mechanics of 
peristaltic transport. Base on this experiment work. 



Shapiro et al. (1971) made a study on peristaltic 
pumping a tube and a channel under the 
assumptions of low Reynolds number and long 
wave lengths. 

Peristaltic transport of chyme modeled by 
a nonNewtonian fluid or power-law fluid in small 
intestine and oesophagus has been studied by 
Srivastava and Srivastava (1985) 

El Shehawey et al. (1999) investigated the 
peristaltic flow of a Newtonian fluid through a 
porous medium. Mekheimer and Al-Arabi (2003) 
have discussed the peristaltic flow of a Newtonian 
fluid through a porous medium in a channel under 
the effect of magnetic field. Elshahed and Haroun 
(2005) discussed peristaltic flow of a Johnson 
Segalman fluid under the effect of a magnetic field. 
Ravi Kumar et al. (2010) studied the peristaltic 
pumping in a finite lengths tube with permeable 
wall. Most of the researchers considered various 
non-Newtonian fluid models under peristaltic 
transport. Ravi Kumar et al. (2011) considered 
power-law fluid in an asymmetric channel with 
permeable walls under peristalsis in their studies. 
Krishna Kumari et al. (2011) studied the peristaltic 
pumping of a Casson fluid in an inclined channel 
under the effect of a magnetic field. 

In view of this, peristaltic pumping of a 
Bingham fluid in a channel with permeable walls is 
studied under long wavelength and low Reynolds 
number assumptions. The effect of various 
parameters on the pumping characteristics is 
discussed through graphs. 

II. MATHEMATICAL FORMULATION 

Consider the peristaltic pumping of a Bingham fluid 
in a channel with permeable walls, under long 
wave length and low Reynolds number 
assumptions (Figure. 1). The flow in a channel is 
governed by Navier - Stokes equations whereas 
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the flow in the permeable wall is described by 
Darcy's law. The channel is of half width 'a '. 



d -dp 



dy 



dx 



(2) 



The region between y = 0 and y = y 0 is called 
plug flow. In this region T ^T 0 . In the region 



between y = y 0 and y = h , we have 



yx 



0 



0 = 



dp 
dy 



ys = 0;ys yy = T n at y = 0 



0 



(3) 



(4) 



The wall deformation is given by 



H (X , t) = a + b sin 



27t 



X 



(X - ct) 



(1) 



where b is the amplitude, X is the wave length 'c' 
is the wave speed. 

Under the assumptions that the channel length is 
an integral multiple of the wave length X and the 
pressure difference across the ends of the tube is 
a constant, the flow becomes steady in the wave 
frame(x,y)moving with velocity 'c' away from the 
fixed frame (X,Y) The transformation between 
these two frames is given by 

x=X-ct, y=Y, u(x,y)=U(X-ct,Y)-c, v(x,y)=V(X-ct,Y) 

where (U, V) are velocity components in the 
laboratory frame and (u, v) are velocity 
components in the wave frame. 

We introduce the following non-dimensional 
quantities: 

_ u x _ y_ h ct ax D 



— 4 1 - pa 2 - Q b — Fa 

lb — — ■ 7) — — ■ q — — (h = — ' F — 



After non-dimensionalisation (after dropping bars), 
the governing equations and boundary conditions 
become 



A \]Da du 
u = ys =-1 —a\y = h 

y a dy 



where y/ is the stream function, Da= 



(5) 



K 



.2 ' 



a 



#=Slip parameter, and T 0 is the yield stress. 

SOLUTION OF THE PROBLEM 

Solving equation (2) and (3) subject to the 
boundary conditions (4) and (5) we obtain the 
velocity as 



2 



LI- 



2 



It, 



h 2 Ph 



-(1- 



h a 

(6) 



We find the upper limit of plug flow region using the 
boundary condition that 



ys y =0 at y = y 0 so we have y 0 = 



Also by using the condition T yx =T h at y — h 



we obtain P = 



_ L h 



h 



Hence 



, — = — = / r,0< / r<l (7) 



h T. 



h 



Taking y = y 0 in equation (6) and using the 

relation (7) we get the velocity in plug flow region 
as 



Ph 2 



^—[(l^+^C^-^-iO^^jo (8) 
2 h a 
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Integrating the equations (6) and (8) and using the 
conditions y/p = 0 at j = 0 , and Iff = If/ at 

y = y 0 

we get stream function as 



Q = -\Qdt = q + \ (13) 



0 



Integrating the equation (12) with respect 
to x over one wave length, we get the pressure rise 
over one cycle of the wave as 



2 / 



Pit . _\_ 

2 ~ Uti 



h 2h a 



(9) 



The volume flux q through each cross - 
section on the wave frame is give by 



yo y 

q = ^Updy + judy 



0 



h 2 (l-f) (10) 



=— [l--(^) + -(^) 3 ]-^ 
3 2 h 2 h a 



The instantaneous volume flow rate Q(x, t) 
in the laboratory frame between the centre line and 
the permeable wall is 



h 



Q(X,t)= \u(X,Y,t)dy 



0 



Ph 



[1- 



2( >0) + 1(^0)3 ] + A 2 

2 h 2 h 



4Da 



a 
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(11) 



From equation (10) we have 



AP = -j 



3(q + h) 



°h 3 [l--T + -T 3 ] + 3h 2 
2 2 



4du 



dx (I 4 ) 



a 



where h(x) = 1 + ^sin 2nx 



The time average flux at zero pressure 

rise is denoted by Q 0 and the pressure rise 

required to produce zero average flow rate is 
denoted by AP 0 so we have 



AP 0 = 



-f 



(3h -3) 



dx 



2 



a 



(15) 



It is observed that when r — > 0 and Da ^ 0 
equation (6), (10) and (14) reduce to the 
corresponding results of Jaffrin and Shapario 
(1971) for the peristaltic transport of a Newtonian 
fluid in a channel. 

The dimensionless friction force F at the 
wall across one wave length is given by 



dP_ 

dx 



-3(q + h) 



2 h 2 h 



a h 



-3{q + h) 



2 2 a 



(12) 
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3(q + h) 
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(l-T) 



Averaging equation (2.11) over one period 
yields the time mean flow Q as 
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against which the pump works. For a given AP , 
the flux Q increases with increasing 1 OC\ For a 



III. DISCUSSION OF THE RESULT 

From equation (14), we have calculated the 
pressure difference as a function of Q for different 

values of x, the ratio of yield stress to the wall 
shearing stress for a fixed amplitude ratio ((/>), 

Darcy number (Da), and slip parameter (OC) is 
shown in figure (2). It is observed that for a 
Bingham fluid, the peristaltic wave passing over 
the channel wall pumps against more pressure rise 
(AP) compared to Newtonian fluid. This type of 
behaviour may be due to the presence of plug flow 
in Bingham fluid. Further, we observe that there is 
no difference in flux for Bingham fluid and 
Newtonian fluids for free pumping case (AP=0). 

For a given AP the flux Q for a Bingham fluid 

depends on yield stress and it increases with 
increasing yield stress. 

The variation of pressure rise with time 
averaged flow rate is calculated from equation (14) 

for different amplitude ratios (<f> ) and is shown in 

figure (3) for fixed T=0.2, Da=0.001, OC=0A. We 
observe that the higher the amplitude ratio, the 
greater the pressure rise against which the pump 

works. For a given AP, the fluxg for Bingham 

fluid depends on yield stress and it increases with 
increasing (j) . 

From equation (14), we have calculated 

the pressure difference as a function of Q for 

different values of Darcy number Da, for fixed 
^=0.6, T=0.2, #=0.1 and is shown in figure (4). 

We observe that for a given AP, the flux Q 

depends on Darcy number and it decreases with 
increasing Darcy number Da. For free pumping the 

flux Q is constant and it is independent of Da. 

The variation of pressure rise with time averaged 
flow rate is calculated from (14) for different values 
of OC and is shown in figure (5) for fixed Da=0.001 , 

T=0.2 and 0=0.6. We observe that the longer 



given flux Q, the pressure difference AP 
increases with increasing OC. 

From equation (16), we have calculated 
the frictional force as a function of Q for fixed 
0=0.6, Da=0.001, OC=0A and for different values 
of T and depicted of figure (6). As T increases 

the flux Q decreases for a given frictional force. 

Further the frictional force for different amplitude 
ratios, different Darcy numbers , and different slip 
parameters is shown in the figures(7)-(9) and it is 
observed that the frictional force F has the 
opposite behavior compared to pressure rise 

(AP). 

We have calculated the pressure rise from 
equation (15) required to produce zero average 

flow rate AP 0 as a function of the amplitude 

ratio 0, Darcy number Da for different values of 

T and is shown in the figure (10). The value of 

AP 0 is larger for Bingham fluid when compared to 

Newtonian fluid (T = 0). As 0— >1, AP 0 
becomes indefinitely large for a given value of T . 
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Plug Flow Region 
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Fig 1 : Physical Model 
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Fig.2: The variation of AP with Q for 
different t with <j> =0.6, Da=0.001, a=0.1 




Fig.3: The variation of AP with Q\or 
different (j) with t=0.2, Da=0.001, a=0.1 
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Fig.4: The variation of AP with Q for 
different Da with ^=0.6, t=0.2, a=0.1 
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Fig.5: The variation of AP with Q for 
different a with ^=0.6, t=0.2, Da=0.001, 




Fig. 6: The variation of F with Q for 
different t with ^=0.6, Da=0.001, a=0.1 
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Fig.7: The variation of F with Q for 
different (f) with t=0.2, Da=0.001, a=0.1 
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Fig. 9: The variation of F with Q for different 
Da with ^=0.6, t=0.2, Da=0.001 
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